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Abstract 

In the framework of the Brueckner-Bethe-Goldstone theory, we determine 
a fully microscopic equation of state for asymmetric and /3-stable nuclear 
matter containing £~ and A hyperons. We use the Paris and the new Argonne 
Av\% two-body nucleon interaction, whereas the nucleon-hyperon interaction 
is described by the Njimegen soft-core model. We stress the role played by 
the three-body nucleon interaction, which produces a strong repulsion at high 
densities. This enhances enormously the hyperon population, and produces a 
strong softening of the equation of state, which turns out almost independent 
on the nucleon-nucleon interaction. We use the new equation of state in order 
to calculate the structure of static neutron stars. We obtain a maximum 
mass configuration with M max = 1.26 (1.22) when the Paris (Avi%) nucleon 
potential is adopted. Central densities are about 10 times normal nuclear 
matter density. Stellar rotations, treated within a perturbative approach, 
increase the value of the limiting mass by about 12%. 



PACS: 26.60.+C, 21.65.+f, 24.10.Cn 
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I. INTRODUCTION 



The nuclear matter equation of state (EOS) is the fundamental input for building models 
of neutron stars (NS) 0. These compact objects, among the densest ones in the universe, are 
produced during the gravitational collapse of massive stars, which explode into supernovae 
at the end of their evolution. Neutron stars are observed as pulsars: because of their fast 
rotation they emit only in particular directions regularly spaced pulses of electromagnetic 
radiation. Although almost 700 pulsars have been detected so far, their gravitational mass 
can be inferred only from observation of a few binary systems ||. The observed NS masses 
are typically » (1— 2)M (where M is the mass of the sun, M & = 1.99x 10 33 g). Above 3M , 
NS are commonly believed to collapse into black holes. Typical radii of NS are thought to 
be of order 10 km, although direct measurements do not exist, whereas the central density 
is a few times normal nuclear matter density (po ~ 0.17 fm~ 3 ). This requires a detailed 
knowledge of the EOS for densities p p . 

This is a very hard task from the theoretical point of view. For instance, the present 
uncertainty on the equation of state at high density implies an uncertainty on the value of the 
maximum mass, important for distinguishing between neutron star or a black hole formation. 
In fact, whereas at densities p ~ po the matter consists mainly of nucleons and leptons, at 
higher densities several species of particles may appear due to the fast rise of the baryon 
chemical potentials with density. Among these new particles are strange baryons, namely, 
the A, £ and S hyperons. Due to its negative charge, the £~ hyperon is the first strange 
baryon expected to appear with increasing density in the reaction n + n — > p + S~, in spite of 
its substantially larger mass compared to the neutral A hyperon (Ms- = 1197 MeV, M\ = 
1116 MeV). Other species in stellar matter may appear, like A isobars along with pion 
and kaon condensations. Moreover, at very high densities, nuclear matter is expected to 
undergo a transition to a quark-gluon plasma ||. However, the exact value of the transition 
density is still unknown because of some technical problems in the QCD lattice calculations 
for finite baryon density. In this paper we disregard these phenomena, because they lie 
outside the scope of Brueckner theory that is applied here. In particular, we concentrate our 
investigation on the production of strange baryons. We assume that a baryonic description 
of nuclear matter holds up to densities as those encountered in the core of neutron stars. 

In a previous article [Q we presented a microscopic investigation within an extended 
Brueckner-Hartree-Fock (BHF) scheme for determining the chemical potentials of the differ- 
ent baryons (n,p, S~,A) in a fully self-consistent manner. We used in our calculations the 
Paris H and the Argonne v 14 nucleon-nucleon (NN) interaction, modified by three-body 
forces (TBF), in order to get the correct saturation point of nuclear matter [0]. In that pa- 
per H we adopted the Nijmegen soft-core M potentials for describing the nucleon- hyperon 
interaction, whereas no hyperon- hyperon potential was taken into account, due to lack of 
restricting experimental data. We mainly concentrated on the calculation of the onset den- 
sity of the E~ and A hyperons. We found that thresholds are reached at densities beginning 
at about 2-3 times normal nuclear matter density, for all the different nuclear equations of 
state considered. 

In this paper we proceed further and present results concerning the equation of state 
of asymmetric and /9-stable nuclear matter containing S~ and A hyperons, obtained in the 
BHF theoretical scheme. We adopt the Paris and the new Argonne v is @ NN potentials, 
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eventually modified by nucleoli TBF according to the Urbana model, and the Njimegen 
nucleon-hyperon potentials. In these calculations still no hyperon-hyperon interaction is 
included. We determine microscopically the chemical potentials of the various species. Their 
concentration is inferred by imposing the conditions of chemical equilibrium, along with 
charge neutrality and baryon number conservation. In general, we observe a softening of the 
equation of state with respect to the pure nucleonic case because of the increased number 
of baryonic species. The main result of our work is that in the presence of hyperons the 
inclusion of the nucleon TBF does not produce any significant change in the equation of state 
with respect to the case with only two-body forces. This is quite astonishing because, in 
the pure nucleon case, the repulsive character of TBF at high density increases the stiffness 
of the EOS, thus changing dramatically the equation of state 0. However, when hyperons 
are included, the presence of TBF among nucleons enhances the population of S~ and A 
because of the increased neutron and proton chemical potentials with respect to the case 
without TBF, thus decreasing the nucleon population. The net result is that the equation 
of state looks very similar to the case without TBF, but the chemical composition of matter 
containing hyperons is very different when TBF are included. In the latter case, the hyperon 
populations are larger than in the case with only two-body forces. This has very important 
consequences for the structure of the neutron stars. Of course, this scenario could partly 
change if hyperon-hyperon interactions were known or if TBF would be included also for 
hyperons, but this is beyond our current knowledge of the strong interaction. 

We apply our new equations of state to the calculation of the static properties of neutron 
stars having hyperonic cores. The values of the maximum mass configuration, i.e., mass, 
radius and central density, turn out almost independent on the chosen two-body nucleon- 
nucleon potential. On the other hand, the values of the central densities are quite different 
when TBF are taken into account, producing less compact stars than in the case without 
TBF. Therefore hyperon stars, i.e., nucleon stars with hyperonic cores, collapse earlier than 
pure nuclear stars. Stellar rotations, which are treated here within a modified Hartle- 



Thorne method [10], increase the value of the maximum mass by about 12%. However only 
configurations rotating at their Kepler frequency have been calculated. 

This paper is organized as follows. In Section we review some formalism in the BHF 
scheme with hyperons, with a discussion of the method chosen for the calculation of the 
baryon chemical potentials. We also discuss the equilibrium compositions among the differ- 



ent baryonic species and the equation of state. In Section |TTT| we illustrate our results. In 



particular, in paragraph [111 A| we discuss the behavior of the single-particle potentials for all 



involved species as calculated in the Brueckner theory. The matter composition is illustrated 
in paragraph I"! [JJ, along with the equation of state. Special emphasis is put on discussing 



the role of the nucleonic TBF in paragraph [III C| . In paragraph [IIID| we analyze the static 
properties of neutron stars and their variations caused by the rotations. Our conclusions 
are finally drawn in Section [TV]. 

II. FORMALISM 

For a detailed account of the extended Brueckner theory including hyperons we refer the 



reader to Refs. [|J and ]LT}|. Here we repeat only the basic formulae and give some details 
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concerning the particular application to neutron star physics. It turns out that in this case, 
besides the nucleons, N = n,p, only the Y = X~,A hyperons appear as stable particles in 
the matter, limiting the baryonic Fermi seas to these four species. (Other types of hyperons 
appear in virtual intermediate states, of course). We remind also that for the present work 
no hyperon-hyperon potentials are used, leading to simplifiction in some places. 

Then, the extended Brueckner scheme requires as input the nonrelativistic nucleon- 
nucleon and nucleon-hyperon potentials. With these potentials, the various G matrices are 
evaluated by solving numerically the Bethe-Goldstone equation |12|], written schematically 



G ab [W]=V ab + J2J2 V *' 



\pp 



Qc 



W - E c + ie 



pp 



G cb [W] 



(1) 



where the indices a, b, c indicate pairs of baryons and the angle-averaged Pauli operator Q 
and energy E determine the propagation of intermediate baryon pairs. In a given nucleon- 
hyperon channels c = (NY) one has, for example, 



E(ny) = mN + my + 
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The hyperon single-particle potentials within the continuous choice are given by 



U Y (k)= J2 U Y N \k)=Re £ E ( kk ' 

N=n,p N =n,p k , <k (N) 

and similar expressions of the form 
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N'=n,p 
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y=s-,A 



apply to the nucleon single-particle potentials. The nucleons feel therefore direct effects of 
the other nucleons as well as of the hyperons in the environment, whereas for the hyperons 
there are only nucleonic contributions, because of the missing hyperon-hyperon potentials. 

These equations ([]]-§]) define the BHF scheme with the continuous choice of the single- 
particle energies. Due to the occurrence of Un and Uy in Eq. @ they constitute a coupled 
system that has to be solved in a self-consistent manner. In our previous work [[| those 
equations were solved for zero hyperon fraction, since we were interested only in their onset 
density. In the present paper we proceed further and perform calculations for arbitrary 
nucleon and hyperon concentrations. Therefore the above equations must be solved for 
several Fermi momenta of the particles involved. 

Once the different single-particle potentials are known, the total nonrelativistic baryonic 
energy density, e, and the total binding energy per baryon, B/A, can be evaluated: 



B 
A 



Pn + P P + Ph- 



Pa 



E 

=n,p,E- 



h Fdkk 2 



Ti 



TXli + 



2rrii 



£nn + £ny 



(5) 
(6) 



4 



with 



E 

N=n,p 



E 

y=s-,A 

E 

y=s-,A 



7T 



d/c fc 2 



7T^ 



A; 2 1 
^ 2mAT 2 



my + 



2m 



7T Z 



^2 

my + h 

Ivtly 



E 

N=n,p 
A") 



k F N) dk k 2 



7T Z 



U?>(k) + U$\k) 



u£-\k) + u£\k) 



(7a) 
(7b) 
(7c) 



Here we have split the energy density into a part due to the action of nucleon-nucleon 
forces, €nn, and due to nucleon-hyperon forces, ejyy. These quantities depend on the total 
baryon density of the system, p = Pn + Py, (pn = Pn + P P , Py = Ps- + Pa), and on 
the baryon fractions Xi = Pi/p, i = p, E~,A. However, due to the fact that the single- 
particle potentials uffi and llffi depend only indirectly and therefore very weakly on the 
hyperon partial densities, the energy density e^N depends to a good approximation only 
on the nucleonic partial density and the proton fraction within the nucleonic subsystem: 
£nn = £nn(pn,% = Pp/Pn), whereas the dependence on the hyperonic partial densities is 
concentrated in e^y. 

This facilitates the determination of the chemical potentials. They are given by 
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where Eq. ([7B) ( J7cD was used in order to arrive at Eqs. ( |8li| , |8"b"| ) (Rc|). 

The contributions to the chemical potentials due to €nn are then the same as in the 
system without hyperons, namely 
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PY=0 



whereas the contributions due to e^y can be expressed by the appropriate components 
of the single-particle potentials, which represent corrections to the Fermi energy of the 
different species. This last step amounts to neglecting certain "rearrangement" contributions 
that appear in the exact expressions for the chemical potentials, namely from Eqs. (0) one 
obtains: 
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and similarly for p and A. Here the last terms represents the rearrangement contributions due 
to the (weak) dependence of on p n and Uy on p^- . In the following, the rearrangement 
contributions of e^y will be neglected, while the much more important e^N is treated exactly, 
as specified in Eq. (H). This simplifies considerably the numerical effort. 

This approximate treatment of e^y is justified, apart from the fact that the nucleon- 
hyperon forces are weaker than the nucleon-nucleon forces, by a peculiarity of the BHF 
approach , where the chemical potential p of a species is given by 



p = e F + U 2 {k F ) + ... . (11) 

Here e F = k F /2m + U(kp) is the Fermi energy, as determined from the BHF single-particle 
potential and U2 is the leading (of second order in the hole line expansion) rearrangement 
contribution to the single-particle potential, which is given by a diagram representing the 
conversion of a hole state into a particle state. It therefore vanishes in pure neutron matter 
for all species different from the neutron, in particular for the proton and the hyperons. 
For the neutron itself, it was shown in Ref. |l3j that in pure neutron matter the second 



order rearrangement contribution is rather small, due to the relatively weak neutron-neutron 
interaction. Therefore, in neutron matter with not too large proton and hyperon fractions, 
the hyperon chemical potentials and the corrections to the nucleon chemical potentials are 
well approximated by the respective Fermi energies. 

As a further simplification we use the fact that in the so-called parabolic approximation 



T4| , the binding energy per baryon in asymmetric (hyperon-free) nuclear matter depends 



to a good approximation quadratically on the asymmetry parameter (3=1 — 2x: 

j(p N , 13) w j(p N , P = 0) + (3 2 E sym (p N ) , (12) 

where the symmetry energy E sym can be expressed in terms of the difference of the energy 
per particle between pure neutron {(3=1) and symmetric {(3=0) matter: 

E sym {PN) = ^{ P n,(3 = 1) - ^{pn,{3 = 0) = 1^^(Pn,(3= 1) • (13) 

The composition of neutron star matter is crucially dependent on the nuclear symmetry 
energy. This quantity strongly affects the onset of hyperon formation, as well as other 
processes like the neutron star cooling rates \TE\ . In the parabolic approximation one obtains 
for the nucleon chemical potentials 

PpApn, (3) = PpApn, 0) - (^(3 2 ±2(3 - I3 2 Pn-^-^ E sym {p N ) , (14) 
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(+ for p, — for n), and in particular 



[fi n - fj, p ](p N ,f3) = 4pE sym (p N ) . (15) 



As far as the hyperon chemical potentials, Eq. flScj), are concerned, in practice an effective 
mass approximation can be employed: 

„ utnr + ° r + £2g!L, (16) 

where the "mean field" U Y = Uy(k = 0) and the global effective mass 

U(k F )-U(0)' 



m 
m 



1 + 



k 2 F /2m 



(17) 



depend on the variables (p, x p , x^-, While U Y depends sizeably on all variables, in 
practice rather good fits to the calculated single-particle spectra were obtained by taking 
into account only the p^ dependence of m*/m. Similarly, the corrections to the nucleon 



Fermi energies appearing in Eqs. ( pa| , pb| ) were found to depend only weakly on the proton 
fraction x, and parametrized as function of pn and xy. 

Once the chemical potentials of all species are known, one can proceed to calculate 
the composition of stellar matter. At high density this is essentially constrained by three 
conditions: i) chemical equilibrium among the different species, ii) charge neutrality, and 
iii) baryon number conservation. The chemical potentials are the fundamental input for 
solving the equations for the chemical equilibrium. At density p ~ po we assume stellar 
matter composed of a mixture of neutrons, protons, electrons, and muons in /3-equilibrium 
[electrons are ultrarelativistic at these densities, p e = (37r 2 px e ) 1//3 ]. In that case the equations 
read 

fin = + , (18) 
He = AV • ( 19 ) 

Since we are looking at neutron stars after neutrinos have escaped, we set the neutrino 

chemical potential equal to zero. Strange baryons appear at density p ~ (2 — 3)po !J> 

mainly in baryonic processes like n + n — > p + S~ and n + n — > n + A. The equilibrium 
conditions for those processes read 

2/i n = p p + p s - , (20) 
p n = p A . (21) 

Further two conditions of charge neutrality and baryon number conservation allow the unique 
solution of a closed system of equations, yielding the equilibrium fractions of the baryon and 
lepton species for each fixed baryon density. They read 

P P = Pe + Pv + Ps- , (22) 
P = Pn + P P + Ps- + Pa ■ (23) 
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Finally, from the knowledge of the equilibrium composition one determines the baryonic 
equation of state, i.e., the relation between baryonic pressure Pg and baryon density p. It 
can be easily obtained from the thermodynamical relation 



od(e/p) 



P B = P'^. (24) 

The total pressure P and the total mass-energy density 8 are then calculated by just adding 
the lepton contributions that are well-known from textbooks, see e.g. Ref. 

P = P B + Pi, (25) 
£ = e + ex . (26) 



III. RESULTS AND DISCUSSION 

A. Single-particle potentials 

In order to illustrate some statements made in the previous section, we show in Fig. [I] a 
representative plot of the single-particle potentials of the different baryons at fixed neutron 
and proton densites, given by pn = 0.4 fm~ 3 and p p / Pn = 0.1, and varying E~ density. 
Under these conditions the S~ single-particle potential is sizeably repulsive, while U\ is still 
attractive (see also Ref. |4j]) and the nucleons are much stronger bound. The S~ single- 
particle potential has a particular shape with an effective mass m*/m slightly larger than 
1, whereas the lambda effective mass is typically about 0.8 and the nucleon effective masses 
are much smaller. 

The influence of increasing S~ density on the hyperonic single-particle potentials is only 
indirect (since there is no hyperon-hyperon interaction) and therefore rather small. There 
is some additional repulsion for the nucleons due to the repulsive effective A£~ interaction, 
growing with E~ partial density. However, these effects represent small variations of the 
single-particle potentials observed in the hyperon-free system, which justifies the approxi- 
mate treatment presented in the previous section. 

The same principal conclusion applies also to Fig. ^, that displays the same information 
as Fig. |l|, but at higher densities p^ = 0.8 fm~ 3 and p p / Pn = 0.2. The main quantitative 
difference is the fact that under these conditions both hyperon single-particle potentials are 
quite repulsive. 

The results displayed were obtained with the Paris nucleon-nucleon potential. With the 
Argonne Vxs one observes only very slight quantitative changes and the corresponding plots 
are not shown. 



B. Stellar matter composition and equation of state 

We come now to the presentation of our results regarding neutron stars. In Fig. |3| we 
show the star composition obtained when the Paris (solid line) or the Argonne V\s (dotted 
line) potential is adopted as nucleon-nucleon force. Three cases are examined, respectively 



8 



(a) no hyperons are present, (b) hyperons are free and (c) hyperons are interacting with 
nucleons. In panel (a) we show the particle fractions versus the baryon density for matter 
containing only nucleons and leptons, ignoring hyperons. We note that neutrons are the 
most dominant species up to very high values of the baryonic density. Electrons and muons 
are present as well. These general features are common to many models. However, the value 
of the proton fraction is model dependent, and it affects strongly the direct Urea cooling 
rates (ITS . 



The picture changes when hyperons are taken into account. This is shown in panels 
(b) and (c). The S~ is the first hyperon to appear, due to its negative charge, whereas A 
formation takes place at higher density. Other hyperon species do not appear in our model. 
In the free hyperon case [panel (b)], the formation of S~ starts at about p ~ 0.4 fm -3 , as 
was published in our previous paper [|]]. In that paper the A onset point was approximate, 
because for its precise determination all the chemical potentials at finite £~ fraction are 
needed. This is performed in the present calculations, and now we estimate exactly the A 
onset point, which is located at about p ~ 0.82 fm~ 3 . As is clearly shown, the thresholds for 
hyperon formation are weakly dependent on the two-body force. The hyperon fractions are 
substantial at high density and constitute a large portion of the stellar core matter. The 
appearance of hyperons induces a deleptonization of the baryonic medium, mainly because of 
the charge neutrality condition. Leptons disappear at high baryonic density, thus hindering 
formation of kaon condensate |L6] . 

When the nucleon-hyperon interaction is taken into account, the scenario described above 
changes quantitatively. In fact, since the E~ Fermi energy is repulsive starting from densities 
just above normal nuclear matter density M, the onset point is shifted to slightly higher 
density. On the contrary, the A formation starts at density lower than in the free case 
because the A Fermi energy is attractive over a wide range of densities. This is clearly 
shown in panel (c). However, the hyperon population is smaller than in the free case and 
deleptonization is less drastic, because the repulsive core of the NY interactions becomes 
relevant at high density. However, even in this case, kaon condensation cannot occur. Again, 
we note how those results are slightly dependent on the two-body interaction. 

The main consequence of introducing additional particle species into matter is the soft- 
ening of the EOS. This softening is essentially due to the conversion of kinetic energy of 
the already present species into masses of the new species. The decrease of the lepton pop- 
ulation does add further softening, but this effect is quite small. The nature of the EOS 
is thus dependent on the number of species, as well as on the details of the strong inter- 
action. It can be instructive to begin with the equation of state corresponding to the case 
(a), i.e., when no hyperons are present. This is displayed in Fig. [| (solid line). On the left 
(right) hand side we show the equation of state obtained when the Paris (Avi$) potential is 
adopted as nucleon-nucleon interaction. That calculation has been improved with respect to 
the one published in Ref. 0], since more channels are now taken into account in the solution 
of the Bethe-Goldstone equation. This has produced a better convergence of the iterative 
procedure for high values of nuclear matter density. 

The properties of the EOS change when hyperons are taken into account. In particular, 
the dotted line corresponds to the case of nuclear matter containing leptons and free hy- 
perons. The presence of hyperons induces a strong softening of the equation of state. The 
inclusion of the nucleon-hyperon interaction produces an equation of state stiffer than in the 
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free hyperon case. This is shown by the dashed line. We observe a similar behavior when 
the Argonne v is potential is used. 

However, it is well known that nonrelativistic calculations, based on purely two-body 
interactions, fail to reproduce the correct saturation point of symmetric nuclear matter 
|jI7| . This deficiency is commonly corrected by introducing three-body forces (TBF) among 
nucleons. This changes the scenario described above. Our method of treating TBF is 
discussed in the following paragraph. 



C. Inclusion of three-body forces 

It is commonly known that a complete theory of three-body forces is not available so far. 
Therefore one has to work with phenomenological approaches. A realistic model for nuclear 
TBF is the so-called Urbana model |L8| , which consists of an attractive term due to two-pion 



exchange with excitation of an intermediate A resonance, and a repulsive phenomenolog- 
ical central term. We introduced the same Urbana three-nucleon model within the BHF 
approach (for more details see Ref. In our approach the TBF is reduced to a density 
dependent two-body force by averaging on the position of the third particle, assuming that 
the probability of having two particles at a given distance is reduced according to the two- 
body correlation function. The corresponding EOS satisfies several requirements, namely 
(i) it reproduces correctly the nuclear matter saturation point 0j, (ii) the incompressibility 
is compatible with values extracted from phenomenology |19[, (iii) the symmetry energy is 
compatible with nuclear phenomenology, (iv) the causality condition is always fulfilled. 

Fig. |5] shows the values of the symmetry energy for the different EOS's that we consider, 
namely the nonrelativistic Brueckner calculations with the Paris and the Argonne vi 8 po- 
tentials with and without three-body forces. For comparison, we report also the symmetry 



energy of a recent calculation performed with a Dirac-Brueckner (DBHF) model [20|, but 
with the Bonn-A potential. We should remind that the DBHF treatment is equivalent to 
introducing in the nonrelativistic BHF the three-body force corresponding to the excitation 
of a nucleon-antinucleon pair, the so-called Z-diagram f22fl , which is repulsive at all densi- 



ties. We therefore expect that the symmetry energy calculated within the DBHF approach 
is always larger than the one obtained in the nonrelativistic BHF calculation, without and 
with three-body force, since in the last case both an attractive and a repulsive component is 
introduced. From Fig. |5| we see that this is true for all the potentials discussed in this paper, 
besides for the Argonne vi$ potential with TBF, maybe because of its strongly repulsive 
core. On the contrary, in the low density region (p < 0.3 fm~ 3 ), both BHF+TBF symmetry 
energies and DBHF calculations are very similar. 

We can proceed now to the discussion of the nuclear matter composition when TBF 
are included in the equation of state. Please note that TBF are included only for nucle- 
ons. Hyperons interact via two-body forces with nucleons (we adopt the Nijmegen soft-core 
potential, as discussed in the introduction), and do not interact at all among themselves. 
Since no experimentally tested hyperon-hyperon potentials are currently available, this as- 
sumption is in line with the exploratory character of this work. Of course, their eventual 
introduction may change the scenarios resulting from our analysis. 

In Fig. ^| we show the matter composition when the TBF among nucleons is included in 
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the BHF calculation. The notation is the same as in Fig. |3|, i.e., solid line represents the 
composition obtained with the Paris+TBF potential, whereas the dotted line corresponds 
to the composition obtained with the Argonne Au 18 potential. Let us first discuss the case 
without hyperons, displayed in panel (a). Because of its repulsive character at high densities, 
the higher value of the symmetry energy allows more easily the conversion of neutrons into 
protons and leptons compared to the case without TBF. The proton fraction can now exceed 
the "critical" value x\j TCSu ~ (H — 15)% needed for the occurrence of direct Urea processes 
0- 

The chemical composition of nuclear matter changes completely when hyperons are taken 
into account. In panel (b) we show the particle fractions obtained when the hyperons 
are free. We notice that the three-body forces shift the onset points of both X~~ and A, 
as previously published 0, down to densities 2-3.5 times normal nuclear matter density. 
Deleptonization takes place, and leptons disappear almost completely just after hyperon 
formation, because now it is more convenient energetically to maintain the charge neutrality 
through S~ formation than /3-decay. At high density nucleons and hyperons are present 
almost in the same percentage. 

This scenario does not change qualitatively when the nucleon-hyperon interaction is 
included, see panel (c). The main difference is that the hyperon onset points move again, 
respectively to higher (lower) density because of the repulsive (attractive) character of the 
S~ (A) Fermi energy at those densities. Even in this case the stellar core is composed by an 
almost equal fraction of nucleons and hyperons. Once again those compositions look very 
much the same both with the Paris and the Avi$ potential. 

The corresponding equations of state obtained when TBF are added to the two-body 
forces are listed in Table | and shown in Fig. |7]. On the left (right) hand side we show our 
results obtained with the Paris (Avi 8 ) potential plus TBF. The solid line shows the equation 
of state of asymmetric beta-stable matter with a percentage of electrons and muons. It 
looks much stiffer than the case when hyperons are introduced. The dotted line represents 
the equation of state of nuclear matter containing free hyperons, whereas the dashed line 
corresponds to the case of hyperons interacting with the nucleonic medium. If we compare 
the equations of state obtained without and with TBF, i.e., the dashed lines of Figs. f| and [7|, 
we see only a small difference. This can be understood by looking at the final compositions, 
respectively panels (c) of Figs. |3| and There we note that the net effect of TBF is a 
decreased presence of neutrons and an enhanced population of A, whereas protons and £~ 
are only slightly affected by TBF, keeping the total percentage of neutral charge almost the 
same in both cases. Therefore, as far as the equation of state is concerned, we expect that 
the relation between pressure and baryon density is very similar without and with TBF. This 
is indeed the case and means that hyperon formation is a mechanism for pressure control, 
as already found by the authors of Ref. |23| within a phenomenological approach. 

Finally, we comment on the deleptonization of the baryonic matter. This is clearly 
indicated by the chemical potentials, drawn in Fig. [8[ In panel (a) we display the chemical 
potentials for neutrons and protons obtained with the Paris potential supplemented by TBF 
(the value of the neutron mass has been subtracted), whereas in panel (b) the electron 
chemical potential is shown. We observe a monotonously increasing function of the baryon 
density when no hyperons are present (dashed lines), whereas the appearance of negatively 
charged hyperons interacting with the medium (solid line) induces a strong deleptonization 
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of the matter. Thus the formation of hyperons and the consequent deleptonization will 

Although kaon and anti-kaon condensation has been 

also recent investigations 



suppress meson condensation |24|,|25 
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proposed as possible state of matter inside neutron stars 
within the relativistic mean field approach p5| find the onset of kaon condensation quite 
unlikely. In any case, the importance of meson condensation will be strongly diminished by 
allowing for the dominant hyperon formation. 



D. Neutron stars 



As already discussed in the introduction, the knowledge of the equation of state is essen- 
tial in order to build models of both static and rotating neutron stars. In fact, the EOS for 



/9-stable matter can be used in the Tolman-Oppenheimer-Volkoff (TOV) equations p6"||2~7 
to compute the neutron star mass and radius as a function of the central density. As a 
first step, we neglect the effects of rotations and calculate the mass-radius relation assuming 
that a neutron star is a spherically symmetric object in hydrostatic equilibrium. Then the 
equilibrium configurations are simply obtained by solving the TOV equations for the total 
pressure P and the enclosed mass m, 



dP(r) Gm(r)£(r) [ + £(Vj 
dr 

dm(r) 



P(r) 



^ i 4nr 3 P(r) 
m(r) 



2Gm(r) 
r 



dr 



47cr 2 £(r) 



(27) 
(28) 



being G the gravitational constant (we assume c = 1). Starting with a central mass density 
£ (r = 0) = £ c , we integrate out until the pressure on the surface equals the one corresponding 
to the density of iron. This gives the stellar radius R and the gravitational mass is then 



M G = m(R)=4n / drr 2 £{r 



R 



(29) 



For the outer part of the neutron star we have used the equations of state by Feynman- 
Metropolis- Teller |28| and Baym-Pethick-Sutherland [29|, and for the middle- density regime 
(0.001 fm -3 < p < 0.08 fm -3 ) we use the results of Negele and Vautherin In the high- 
density part (p > 0.08 fm -3 ) we use alternatively the equations of state displayed in Table |. 
For comparison, we also perform calculations of neutron star structure for the case of purely 
nucleonic asymmetric and /3-stable matter with some lepton fraction. 

The results are plotted in Fig. |9|, where we display the gravitational mass M G (in units 
of the solar mass M & ) as a function of the radius R and the central baryon density p c . 
We report the results obtained using either the Av is or the Paris two-body potential with 
three-body forces. The solid lines represent the equilibrium configurations of neutron stars 
composed only of nucleons and leptons; the dashed lines show the configurations of stars 
whose composition includes interacting hyperons. As we can see, the softening of the equa- 
tion of state due to the presence of additional baryonic species produces a strong decrease in 
the value of the limiting mass, from about 2M down to 1.2M , almost independent on the 
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nucleon-nucleon interaction. The limiting central densities stay nearly constant at about 7 
times normal nuclear matter density. 

For clarity we report in Table [TI| the properties of the maximum mass configuration 
obtained with equations of state of /3-stable asymmetric nuclear matter with and without 
hyperons. Without hyperons we observe mainly the effect of three-body forces among nucle- 
ons on the equation of state and, therefore, on the values of the maximum mass configuration. 
In particular, the increased repulsion among nucleons produces a stiffer equation of state 
and a higher value of the maximum mass, with a smaller radius and central density. Even 
in this case we note a strong similarity in the limiting values independently on the two-body 
potential. 

The situation changes dramatically when interacting hyperons are included. In this case 
the additional repulsion produced by the nucleon three-body forces is counterbalanced by 
the increased population of hyperons, thus leading to a very soft equation of state. The 
values of the maximum mass do not differ very much from those calculated without TBF. 
However, the values of the maximum central densities are substantially different, despite 
the EOS looking similar in the two cases. When TBF are included, this value becomes 
smaller, see Table ||. Without TBF, in fact, the EOS is slightly stiffer at higher densities 
and the corresponding curve of the mass vs. central density (not reported in Fig. 9) has a 
shallow maximum at larger central density. However, the corresponding maximum mass is 
only slightly larger than in the case with TBF, as reported in Table [H]. This also shows the 
sensitivity of the results on the details of the nuclear EOS. It has to be noticed that, for 
stable configurations, at a given value of the neutron star mass the central density is larger 
when TBF are included. Furthermore, neutron stars built with equations of state of matter 
containing baryons interacting via TBF possess a core with a larger hyperon population 
than in the case without TBF. 

The above scenario changes when the rotations are included. In particular, in order to 
treat stellar rotations in general relativity, we follow the method discussed in Ref. |]10|| . In 
those papers the Einstein equations for rotating massive objects are solved perturbatively 
with a modified version of Hartle's method. With this method the authors are able to 
determine sequences of star models rotating at their respective general relativistic Kepler 
frequency Qk- in Fig- we display the net effect of rotations on the mass-density relation. 
There we display star sequences at equilibrium both for static (dashed lines) and rotating 
(solid lines) neutron stars. The upper curves refer to equations of state without hyperons, 
whereas the lower curves refer to equations of state containing hyperons. The rotations 
increase the value of the limiting mass by about 12%, decreasing the limiting value of 
the central density. Again, we do not observe a significant difference for the various NN 
potentials. In Table |H] we report also the values of the limiting configuration. 

IV. CONCLUSIONS 

The principal finding of this paper has been a suprisingly low value for the maximum 
mass of a "neutron star" that barely comprises the "canonical" value M ss 1.4M . While 
certainly the technical insufliciences and approximations of our approach can account for an 
uncertainty on this limit of a few percent, it seems rather difficult theoretically to avoid a 
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rather low limit on the maximum mass. As we have seen, the mere presence of additional 
baryonic degrees of freedom in the form of hyperons renders the maximum mass quite 
insensitive to the stiffness of the nucleonic equation of state. Because the hyperonic onset 
densities at about 2-3 times normal density seem to be rather robust and model independent, 
it seems that the maximum mass can be substantially increased only if the nucleon-hyperon 
and/or hyperon-hyperon effective interactions become extremely repulsive at high density. 
Unfortunately at present the experimental information on the nucleon-hyperon interaction 
is rather scarce, and for the hyperon-hyperon case practically non-existent. 

Nevertheless, recently the Nijmegen NY potential model was extended to the full YY 



case by imposing isospin symmetry ||31|| . In Ref. |32[ these new potentials have been used to 
include the hyperon-hyperon interaction in neutron star studies within a similar approach, 
using a variational nucleonic EOS together with hyperon single-particle potentials deter- 
mined in a standard choice Brueckner scheme. The obtained maximum neutron star masses 
turn out to be quite similar to the ones presented in this work. 

Let us also mention that, within different versions of the relativistic mean field theory, 
the introduction of hyperons in the EOS of nuclear matter leads to maximum masses that 
are only slightly larger [|33| or substantially larger [34j than the ones reported here. This 
illustrates the lack of constraints imposed on this type of models. 

In Ref. [33 1 the possible onset of quark-gluon plasma in the interior of neutron stars is 
also investigated and found to be unlikely to occur up to baryon densities we found in the 
interior of the neutron star. It is expected that the possible presence of the deconfined phase 
would produce an additional softening of the EOS, and therefore a further decrease of the 
critical mass. 

It seems, therefore, that indeed the theoretical limit of the maximum mass of a neutron 
star is rather small. A lowering of the maximum mass of neutron stars due to the onset 
of hyperon formation was actually already qualitatively anticipated in the very first articles 
which investigated baryonic matter with hyperons J55|. However, no quantitative analysis 
could be done at that time, due to the uncertainties in the theoretical models. Within the 
present microscopic approach we can give now a more quantitative prediction of this effect. 
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TABLES 



TABLE I. EOS for /^-stable matter with hyperons obtained in the BHF approximation using 
the Argonne v±s or Paris two-body interaction complemented by the Urbana model for three-body 
forces. We display the baryon density p, the total mass density £, and the total pressure P. 



o ( fm 


£ no 14 1 


Paris+TBF 


P (10 34 
Au-is+TBF 


flvn cm ^\ 

Paris+TBF 


0.08 


1.27 


1.27 


0.07 


0.06 


0.16 


2.64 


2.63 


0.42 


0.39 


0.2 


3.33 


3.33 


0.80 


0.75 


0.3 


5.13 


5.12 


2.62 


2.52 


0.4 


7.11 


7.20 


4.70 


5.05 


0.5 


9.01 


8.96 


7.16 


7.40 


0.6 


11.1 


11.0 


9.71 


10.1 


0.7 


13.0 


13.1 


12.2 


12.8 


0.8 


15.3 


15.2 


15.3 


16.2 


0.9 


17.3 


17.2 


18.3 


18.9 


1.0 


20.0 


20.1 


22.5 


23.5 


1.1 


21.8 


22.1 


25.7 


27.2 


1.2 


23.9 


24.4 


29.3 


30.5 


1.3 


26.0 


26.2 


33.3 


34.5 


1.4 


28.8 


28.5 


38.8 


39.2 


1.5 


31.1 


31.1 


43.7 


45.2 


1.6 


33.5 


33.4 


49.0 


49.7 


1.7 


35.4 


36.2 


53.2 


56.2 



TABLE II. Properties of the maximum mass configuration obtained for different equations of 
state: Mq is the gravitational (maximum) mass, R is the corresponding radius, and p c the central 
baryon density. In each case the results of the EOS without hyperons (no Y), including hyperons 
(Y), and including rotation at the Kepler frequency (Y+Rot) are listed. 



EOS 




M G /M Q 






R (km) 






Pc (fm 3 ) 






no Y 


Y 


Y+Rot 


no Y 


Y 


Y+Rot 


no Y 


Y 
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1.64 


1.26 


1.44 
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8.70 


10.53 


1.53 


1.86 


1.56 


Paris 


1.67 


1.31 


1.50 


8.90 


8.62 


10.30 


1.59 


1.84 


1.62 


Av 18 +TBF 


2.00 


1.22 


1.41 


10.54 


10.46 


12.68 


1.11 


1.25 


1.09 


Paris+TBF 


2.06 


1.26 


1.45 


10.50 


10.46 


12.65 


1.10 


1.25 


1.09 
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FIGURES 



p N = 0.4 fm 3 , p p /p N = 0.1 , Paris NN & Nijmegen NY 
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FIG. 1. The single-particle potentials of nucleons n, p and hyperons £~, A in baryonic matter 
of fixed nucleonic density p^ = 0.4 fm -3 , proton density p p /pn = 0.1, and varying £~ density 
Pe- / ' Pn = 0.0,0.1,0.2,0.3. The vertical lines represent the corresponding Fermi momenta of n, p, 
and For the nucleonic curves, the thick lines represent the complete single-particle potentials 

— (n) (p) 

Un, whereas the thin lines show the values excluding the S contribution, i.e., + Upj . 
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FIG. 2. Same as Fig. 0, but for p N = 0.8 fm" 3 and p p /p N = 0.2. 
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FIG. 3. Composition of neutron star matter in the BHF model. The top panel (a) shows the 
results without hyperons; the middle panel (b) results with noninteracting hyperons; and the lower 
panel (c) includes interaction between nucleons and hyperons. 
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FIG. 4. EOS of neutron star matter in the BHF model. Results without hyperons (solid 
lines), with free hyperons (dotted lines), and with interacting hyperons (dashed lines) are shown. 
The Paris (left) or Argonne vis (right) nucleon-nucleon potentials were used. 



19 





20 




Baryon density p (fm 3 ) 



FIG. 7. As Fig. ||, but including three-body forces. 
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FIG. 8. Chemical potentials of the different species as a function of baryon density in neutron 
star matter. The upper panel (a) shows results for neutrons and protons, whereas the electron 
chemical potential is displayed in panel (b). The dashed (solid) line denotes the case of matter 
composition without (with) interacting hyperons. 
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FIG. 9. Neutron star mass as a function of radius (left) or central baryon density (right) in 
the BHF+TBF model. Results without (solid lines) and with interacting hyperons (dashed lines) 
are compared. 
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FIG. 10. Neutron star mass as a function of central baryon density in the BHF+TBF model 
with and without hyperons. Curves for nonrotating stars (dashed lines) and rotations at the Kepler 
frequency (solid lines) are shown. 
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